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ABSTRACT 

^ , It is shown that for a large class of non-holonomic quantum mechanical systems one 

QQ ■ can make the computation of BRST charge fully algorithmic. Two computer algebra 

|/~\ ' programs written in the language of REDUCE are described. They are able to realize 

f*"^ the complex calculations needed to determine the charge for general nonlinear algebras. 

Some interesting specific solutions are discussed. 
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1. Introduction 



Q^l BRST theory^ has proven to be a very powerful tool to treat gauge theories and, 



in particular, to quantize gauge field theories. Indeed it may be useful to get the 
physical states from a space larger than seems necessary, because the restriction of 
the phase space to the physically meaningful configurations often hides the natural 
symmetries. Some systems, such as Yang-Mills theory, do not even admit global 



> 

X, 

H , gauge fixing conditions. In such cases BRST theory naturally introduces the ghosts, 

necessary to write the path integral formulation of gauge theories. 

At the core of BRST theory is the BRST charge Qb, used to select the physical 
states from a larger space, and specify the physical equivalence of two a priori 
different states. 

QB\'^phys>^0, (1.1) 

l^phys > « l*;,,^;, > +Qs|*0 > • (1.2) 
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The BRST charge is closely related to the symmetries of the system described 
in the Lagrangian formulation through Noether theorem. Its construction can, 
however, be based only on the knowledge of the algebra of constraints. 

A recursive procedure to build the BRST operator order by order in the ghosts 
has been developped but that procedure is not fully algorithmic. We show here 
that the procedure can be made algorithmic for an important class of theories. We 
describe the capabilities of two programs written in REDUCE^ which do indeed allow 
to compute the corresponding BRST charge. This class covers quadratic algebras 
such as those discussed in ref. 3,4, thus extending the realm of application of BRST 
theory beyond usual gauge theories. The use of computer algebra is mandatory for 
all but the most simple cases as the computations involved are tedious, although 
systematic. 

The paper is organized as follows. 
In sec. 2, the construction of the BRST operator is briefly described and two algo- 
rithms for its construction are given. The programs are described in sec. 3 while 
their applications are illustrated in sec. 4. Conclusions as well as the present limi- 
tations of these programs are presented in sec. 5. 

2. The BRST operator 

In the Hamiltonian formulation of a gauge theory, the symmetries of the system 
arc described by the constraints, which are relations between the generalized coor- 
dinates and the momenta. The constraints of the Hamiltonian formulation are not 
all related to gauge invariance. Those related to it are called first class constraints. 
The others, said to be of second class, can be eliminated by an adequate redefinition 
of the brackets. 

Let us consider a set of N first class constraints 

{Ca, a^l,...,N}. 
They satisfy a Poisson bracket algebra 

{Ca,Ct}^L,'C, (2.1) 

where the f's, the structure functions depend on the canonical variables. Here and 
in the following, summation over repeated indices is understood except otherwise 
stated. 

The computation of the BRST charge can be made fully algorithmic if we assume 
that one can write the structure functions in terms of the constraints. Eq (2.1) can 
then be written as 

{Ca,a}^j2f^a,^^"I[^^. (2-2) 

where, now, the /^'s are constants and g is a number obtained from the relations 
between the constraints and the structure functions. In the present work, the nec- 
essary manipulations to obtain it are not considered so that Eq. (2.2) is the starting 
point. 



To each constraint Ca , one associates two variables of odd Grassniann parity : 
a ghost ?7° with ghost number 1 and a ghost momentum Va with ghost number -1. 
They satisfy the following bracket algebra"^. 

K,^'} = {7'a,n} = o, K,n} = -<5?. (2.3) 

The BRST operator fi is defined by 

- its nilpotency ({J7, 51} = 0), 

- its ghost number 1 

- it involves the term rj'^Ca 

It can be decomposed in powers of the ghosts as 

n^Y. "^"^ = E K.:i:^y' ■ ■ ■ v'-^'Tn ■■■'Pj. (2.4) 

n=0 n=0 

where the ^'s depend only on the constraints and structure constants of the prob- 
lem. It can be determined in a recursive way by 

n^^^^V^'Ca, (2.5) 

Sfl(n+1) = _£)(«) (2.6) 



with 



^F = -F^a (2.7) 



and 

n-l 



\fc=0 fe=0 / 

The second term is not present when n = 0. Here, {•, ■}c and {•, ■}gfi are respectively 
the constraint and ghost Poisson brackets. 
The computable D is written in the form 

jjin-l) ^ 2ia)a,-a„^,b^,^^^^^^ . . . p^^_^Cb (2.9) 

in order to make the computation algorithmic. The Z's depend on the ghosts 77"^ 
and structure constants and are antisymmetric over the ai, • • • , a„_2, b. (a) is a set 
of indices which may be empty. Note that this expression is not easy to obtain 
explicitly, as it involves a decomposition of a polynomial (D^""^') in products of 
the constraints, seen as generators of a polynomial ideal. 
From Eq. (2.4), one writes 

r!(") =if(")"i-""-iC(„)7'ai •••K_i (2.10) 



"The choice of sign in the bracket of ghosts and ghost momenta is conventional. In this paper, 
the conventions of ref. 1 are adopted. 



where x^'^^"^'"'^"-^ is antisymmetric over ai, • • • , a„_i. C(^a) denotes the product of 
constraints whose indices are in the set (a). 
Eq. (2.6) leads to 

(n - l)x(")''--'^'-^''C(„)7'a, • • •T'a^.Cfc = z(")"--""-^^C(,)P,, • • -Pa^.Cb. 

(2.11) 
Its solution is 

j^(a}ai---a„_2b _ ( ^(a)af-an-2b i g(a)ai ■■■a„_2b ) (2 \2^ 

n — 1 \ / 

as Z is antisymmetric as mentionned. 

The tensor G is arbitrary except that it is antisymmetric over the set of indices 
formed by the union of ai, • • • , a„_2, h and one index of the set (a). In practice, it 
will be built from the tensor Z with the above antisymmetrisation, each possible 
term being multiplied by an arbitrary coefficient. 

Since the BRST charge is defined only up to an arbitrary BRST exact term, 
there exist infinitely many solutions of Eq. (2.6). In particular, the Q,^^> can be 
built as an arbitrary linear combination of the type 

n-l 

^^^^ai...a._i(a)a,...a„_i ^2.13) 

i=0 

where no antisymmetrisation on the indices of Z is involved. Some of the cocfhcicnts 
ai will be fixed by imposing Eq.(2.6). It is left to find heuristic criteria to fix the 
remaining ones. Two of them are explained in the next section together with the 
corresponding programs. 

3. Description of the Programs 

A preliminary step to the writing of the programs is to find ways to fix all 
constants a^'s. This is done in two ways. They are successively described. 

In the first one, during the calculation of each $7^"^ the coefficients a^'s are 
chosen to minimize the number of terms. 

In the second one, one assumes that the BRST charge remains invariant under 
a redefinition of the constraints, of the ghosts and the ghost momenta. This kind of 
constraint is usual in homological perturbation theory. The principles at the basis 
of this choice are the following: 

Suppose one can define a linear operator cr, and a partition A — (BkAk of the 
algebra of the constraints such that, for all n, 

{Sa + aS)A = NnA (3.1) 

for all A e An- Such an operator is called a contracting homotopy operator for N. 



Now consider the equation (2.6) where £)(") is known, and f7("+^) has to be 
determined. Decompose them as 

^(„+i) ^ j2 n["+^^ f]("+^^ e A, (3.3) 

i 

i?(")=^i^(") i^(")eA. (3.4) 

i 

5Di'' — since SD^'^'' = and 6 is hnear. Thus, solving Eq. (2.6) is equivalent 
to solving Sriy^ ' ~ —Di for all i, as 

2 2 i 

(3.5) 
Now, since 

JaiPf ^ = ((5a + a(5)i?f ) = N^D^^''^ , (3.6) 

as (5D,["^ = 0, one sees that 

^^"+'^ = ^^^r' (3.7) 

is effectively a solution at degree i. 

One can then state that for a contracting homotopy a for N, 

"^"^'' = E ]J--^1"' (3-8) 

is a solution of the equation (2.6). 

A particular contracting operator a is given by 

. = ^r. (3.9) 

where the derivative with respect to the constraints is defined since the expressions 
considered depend only on the phase space variables through the constraints. 

This definition of the contracting operator corresponds to a particular choice of 
the values of the parameters a^ . Other choices would result in other partitions of 
the set of expressions considered, and other contracting operators. 

The BRST charge resulting from the current choice of contracting homotopy 
operator is invariant under linear redefinitions of the constraints: 

Ca^C', = Aja (3.10) 

with the appropriate redefinitions of the ghosts and ghost momenta: 

Va -^ V'a = Aj'Vb 

f^'^ ^ ry'" = (A-l)b«r,^ (3.11) 



One is ready now to describe the two programs computing the BRST charge. 
They are based on the two choices of constants ai 's described above. On the other 
hand, one is written in the algebraic mode of REDUCE, the other is written in its 
symbohc mode. A subsidiary usefulness of writing two programs is to control the 
validity of the calculations. Apart from the algorithms, their originality does not 
come from the handling of anticommuting variables (the ghost and ghost momenta) 
but from the extended use of dummy indices to represent the various expressions. 
Thanks to that, one is able to apply them without reference to the number of 
constraints and, of course, the complexity of the computations is also independent 
on this number. Expressing the brackets {ry^Ca, AT} for example involves as many 
calculations as there are constraints if the summation is explicit. On the other hand, 
non explicit summations on dummy variables allow the treatment of generic cases. 
However, the problem is now to achieve full simplification of polynomials. Outside 
the context of tensor algebra, there exist no package which can do that except the 
package DUMMY^ recently created by one of the authors. Both programs also use the 
package ASSIST^ but they use different functionalities included in it. 

3.1. Program I 

Input: 

- Algebra of constraints Eq. (2.2) 

- rj(") Eq. (2.5) 
Output: f7(") 

The process of computations is: 
fc = 
LOOP: 

compute Z?*^*^) from Eq. (2.8) 

extract the corresponding Z from Eq. (2.9) 

construct rj^'^+i) from Eqs. (2.12) and (2.10) 

arbitrary coefficients left in solution optionally fixed as explained below 

if A: = n - 1 then return fi(") else k:=:k+l go to LOOP 

In the program one has to handle polynomials of the type 

Y^ati::Z^-''-Ca,...Ca,v'''---v''-'Pc,...Vc^ (3.12) 

where 77 and V are odd Grassmann variables and the a's stand for Z or K . It is useful 
to consider the ghost and ghost momenta products as antisymmetric operators : 

r;''^.../;''™ -77^i-^™, 7'ci...Pc„-7'ci...c„. (3.13) 

a. calculation of D 

Brackets of expressions like (3.12) can be decomposed into two parts, the con- 
straint brackets 

\ '\ ' „ai---a;Ci...c„iai---",/Ci...c^, J- n n r<' \ 



^6i...6„6',..Op^^ ^^^, ^^^(_1)"™' (3.14) 

and the ghost brackets 

\ ^\ ^ „ai---a;Ci...c„iai---",/C^...c^, r< n r< 

2^ 2^°'bi...b^ b[...b'^, ^ai ■ ■ ■'~'a,^a[ ■ ■ ■'-'a,, 

{r/''-''"7'c,...c„, /-■■'''"' ^ci..<,}. (3.15) 

In both cases, the time of calculation is considerably shortened if either the C's 
or the ghosts or both factorize. In practice, a factorization of the ghosts is more 
frequent than a factorization of the constraints. 

Using the Leibniz rule, the constraint brackets are easily computed from the 
input. 

The ghost brackets implement the formula 

\^\^ ai...ak-iak + i...ai-^a[...a'-^ /'_ 1 ^fe+^+l+i2 jl X^fc , 

Z^Z^n I bi...bi^b[...b[__^b[^^...b'.^y i-) Ofoj + 

.fc=i 1=1 

31 12 

V^ V^ ai---a>ia'i.--aL-l"'fc + i---"ji-Pi ^ ^ i i, i, K'^'k ( _-\\k+l+3lJ2+ili2+i2Jl 

2_^2^n ^ I bi...bi_ibi+i..M2b{..M^"bi \ '-) 

k=l 1=1 

(3.16) 

This formula can be checked by any program dealing with anticommutating vari- 
ables. The D function is obtained from the bracket calculation. Because it involves 
many terms with the same structure, the use of the function CANONICAL from the 
package DUMMY is essential to simplify it. Because (JD^") = 0, one checks this vanish- 
ing for each value of n. It is interesting to note that, in all the considered examples, 
5D^^' — gives all the Jacobi identities for the structure constants of the algebra. 

b. determination of the Z function 

Because the result for D is obtained from CANONICAL, an antisymmetrization 
over the ghost indices is necessary to extract Z. After this antisymmetrization, 
each monomoial is divided by the constraints and the ghost momentum operator. 
This gives Z. 

c. construction of the BRST operator 

It is constructed from Eq. (2.12) where the arbitrary G"s are built from Z by 
a further antisymmetrization and multiplication by arbitrary coefficients. These 
appear in the output but can also be fixed if one requires the number of terms 
appearing in the output to be minimal. 



3.2. Program II 

The input and output for the second program are identical to that of the first 
program: Input: 

- Algebra of constraints Eq. (2.2) 

- rj(") Eq. (2.5) 
Output: f7(") 

The process of computation is also similar, as it follows the standard steps 
described earlier. 

The BRST charge computed here differs from that returned by the first program 
in the choice of values for the coefficients a. We have indeed chosen to adopt here 
the contracting operator (3.9) which yields a BRST charge invariant under linear 
transformations of the constraints. 

This choice induces a larger number of terms in the expression of 17, and thus re- 
quires more computer resources. It was therefore necessary to work in the symbolic 
mode of REDUCE. 

A further particularity of this program compared to the previous one lies in its 
handling of the Jacobi identity. Indeed, these were useful in reducing as much as 
possible the number of terms in the expressions without changing the contracting 
homotopy operator. We have built a procedure returning a normal form of polyno- 
mials with respect to the Jacobi identity for the particular algebras studied. This 
enabled us not only to reduce the size of the expressions, but also to check the 
validity of the results returned. 

It should be noted however that this handling of side relations is very time 
consuming, and requires ad hoc tweaking for each new algebra. 

Finally, this program has been used to study a partial classification of polynomial 
Poisson structures^. 

4. Results 

Various algebras have been considered 
Usual linear Lie algebras 

The constraint algebra is given by 

Use of both programs gives 

^^'^ - \fabV'''V.- (4.2) 

Computing D^^^ and checking the SD'^^^ = lead to the Jacobi identity 

/[a6/|c-0. (4.3) 



This leads to Z?^^^ = and stops the construction because aU the added contribu- 
tions vanish. 

Self-reproducing algebras 

Let us consider the algebra 

{Ca.Cb}^Tal>CaCb (4.4) 

where no summation over a,b is involved. In this section there is no summation over 
repeated indices. Such an algebra is characterized by the fact that Jacobi identities 
are trivial. At each order from the second one, an arbitrary coefficient is generated. 
The number of handled terms increase quickly. Fortunately, a particular choice of 
the arbitrary coefBcient a„ occuring at the order n by 

a„ = (4.5) 

n 

considerably simplifies the result and the BRST operator can be obtained with 
program I in a closed form. 

"-E IJ T.iK»rc>,v' (4.6) 

n=0 ■ b 

where 

ifh = ^?7"^ar,6 (4.7) 

a 

The program is also able to compute 51 for arbitrary coefficients. An example of 
output is given in Appendix A. The expressions are much more complicated but 
they allow to check the compatibility of the calculations made by both programs. 

Pure quadratic algebras 

The BRST operator has been computed for the constraint algebra 

{Ca,a}^dltCcCd (4.8) 

along the same lines as above. One gets 

^^'^ = yHv^'VcC,. (4.9) 

The vanishing of SD^^'' gives the Jacobi identity 



<.^<?-0. (4.10) 



The calculation of fl^^^ gives 



"^'^ - ldli4^v'"''PafC9- (4.11) 



Program II allows to easily compute Q up to order six. 

Mixed linear and quadratic algebras 

The following algebra 

{Ca,Ch} = U^C, + dliC,Cd (4.12) 

with 

dltdfj = (4.13) 

is an extension of an algebra studied by Schoutens, Sevrin and Van Niewenhuizen^. 
Applying the programs, one gets 

^^'^ = I {.fab'v'"'r. + ditv^'r^c,) . (4.14) 

The vanishing of SD^^'' again leads to the Jacobi identities 

flJfd]c - 0, (4.15) 

fla^dll + df^J^J + dllj,f^ = 0. (4.16) 

Together with (4.13), they imply D^^) = and f^^^) ^ q. 
Jl(3) ig given by 

^^'^ = ^dZdfj^^^r'^'r.f,. (4.17) 

One can note that here and in the previous example, no arbirary coefficient 
is involved. The condition (4.13) implies the vanishing of higher orders. If this 
condition is released one can again compute easily Q up to order six with program II. 

An example of cubic algebra 

The cubic algebra generated by the C's 

{C,,,CX} = .f^:,,C,, + Dif,lC,,C,, + E'/JfCd,C,,C,, (4.18) 

is an extension of the wcllknown spin 4 algebra. The various structure constants 
f,D,E are antisymmetric over their lower indices and symmetric over their upper 
indices. They satisfy the Jacobi identities : 

flJfd]c = 0, (4.19) 

.hJdli + DI:J^/^ + Dlij,^l ^ 0, (4.20) 

2DJ,'^I^^ D'^/Y; + EJ.'^l^"^'^ f^-', , + S&' E'^"^'}'^' = 0, (4.21) 

^j^{d,d, ^d,d<,}d^ 2El,'^f;l"^'- D'^fY' = 0, (4.22) 

d-i_[d2 djds,] di[d2 d-jdg,] ' ^ ' 

Et^'-Ettt^ = 0. (4.23) 



The BRST charge for this algebra is given to order six in Appendix B. 

5. Conclusions. 

As shown by the results in the previous sections, the new programs allow the 
computation of the BRST operator when the algebras of constraints are more com- 
plicated than the usual linear algebras. We recall that quadratic algebras have been 
considered recently in the framework of the study of superconformal field theories'* 
Program I is written in the algebraic mode of REDUCE . It is, of course, less efficient 
than the program written in symbolic mode but is still quite able to make the most 
relevant computations in a reasonable time. Its main limitation is the fact that it 
does not take properly into account the Jacobi identities. For instance, the calcula- 
tion of the first nontrivial term after n^^' is incorrect by a numerical factor g. The 
necessity to take into account the Jacobi identities can be implemented in program 
I but program II fulfills that job and, so it does not look to be worth the task. As 
far as this aspect is concerned, it should be stressed that it is because of the high 
level of symmetry of Poisson algebra structures that this was possible and it is not 
claimed that the algorithm is efhcient. Finally, the specificity of program I with re- 
spect to program II lies in the calculation of the BRST operator in self-reproducing 
algebras. The result can indeed be written in a more compact form than the choice 
of the contracting homotopy made in program II allows to. Programs are available, 
on request, by electronic mail at the following address caprasse@vml.ulg.ac.be. 
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Appendix A 

omega{3) computed from program I for arbitrary parameters alpha2, alphaS : 
omega{3) :— brstconstr{3,any); 

omega{3) := (alop{sl, s2) * alop{sl, s3) * eta(sl, s2, s3, s4)* 

(—3 * alop(s\, s4) * contr{sl) * prond{s2, s3, s4) * alpha2 * alphai 

— 3 * aZop(sl, s4) * contr{sl) * prond{s2, s3, s4) * alphaZ 

+ 9 * alop{s\, sA) * contr{s4) * prond{sl, s2, s3) * alpha2 * alpha3 

+ 12 * alop{sl, s4) * contr(s4) * prond{sl, s2, s3) * alpha2 

+ 9 * aZop(sl, s4) * confr(s4) =¥ prond{sl, s2, s3) * alphai 

+ 12 * alop{sl, s4) * contr(s4) * prond{sl, s2, s3) 

+ 56 * alop{s3, s4) * contr(s3) * prond{sl, s2, s4) * alpha2 

+ 84 * alop{s3, s4) * contr(s3) * prond{sl, s2, s4) 

+ 8 * alop{s3, s4) * conir(s4) * prond{sl, s2, s3) * alpha2 

+ 12 * alop{s3, s4) * contr(s4) * prond{sl, s2, s3)))/96$ 

omega[3) computed from program I with alpha2 = —3/2, alphaS = —4/3 : 
omega{3) :— brstconstr {3, simplify); 

omega(3) := {alop(sl, s2) * alop{sl, s3) * alop{sl, s4) * alop{s\, s5) * alop{sl, s6)* 
conir(sl) * eto(sl, s2, s3, s4, s5, s6) * prond{s2, s3, s4, s5, s6))/3840 



Appendix B 

An example of the output from Program II for the cubic algebra follows. The formula 
below is a direct TJiJX output of package TRI of REDUCE. It should be clear that we can 
compute n to higher orders. 

f^(6) = (ry'^i ry'^^ r^"^ rf^ ry<^ rf- rf^ Vd, Pd, Vd^, Pd,, Pd,, Pd,, 

°^^Jdi4d2 -Idisdie J dndis ^<iia ^d^di ^d^de 

j^di4disds I Qf^A fdg fdio rdii y-1 T^di7digdi2 

^djdi ^ °"^.'di4di5 Jdi6d2 J dndis ^^19 ^d^di 

Tpdiidigdis Tjjdi^digdg \2(} f^^ fdio r)di8di2 

d^de djd\ 'I di^d-i^ ■> diQdir d^d^ 

-TidiQdii f-i rpdi4dirdis rpdi^digds fio {dg r^di^dio 

^dwd2 '-'dia^dsda ^djdi ^'^ Jdiidis-^d2d3 

j~)dndii j~)di8di2 ptdudia fi rpdi^digds _ tn rdg 
^did5 ^disdn-^disde '-^dig -t^drdi ^^ J d^dis 

-jdiedio 



^d2d3 ^did^ ^dieds ^disdn '^d^g J^dydi ^"° 

fdg T~)diedio ]~)dndii T~)di8di2 T~)di4,di3 (~i rpdi^digdg 

Jdiiditi ^d2d3 ^diedi ^d^d^ ^ d^ds '-'dig ^djdi 

I QA r)di4dii j-)dijdi2 r^tZigrfis p^di^dio p^digdg r^di^dg f^ 
~^°^^did2 ^didi ^dsde ^ dudi^ -^ dndi^ ^ digdr ^die 

I FC*^ r^di4dii j-^di7di2 j-^digdi3 T-jdigdio f^di^dg j-^disds /-i 
^"^"-^^1^2 ^d3d4 ^dr-,de ^diidi^-^dndm digdT^dia 

I 1 ao r^dijdii r^disdi2 pidigdi3 T-jrfie^io pidi^dg T^di^^d^ f^ 
+ ^"°-^dirf2 ^d3di ^dsde -^di4di6 -^diTdis di9(i7 ^"^le 
-i- 99/1 T^di4dio j-^dirdii p^di^dg T-jdigdg j-\digdi3 T~)rfi5di2 /^ 

-r ZZ'iiy^^^^ -^^4^5 ^di4di5 -^diydi -^dis^T ^digdg ^di6 

I O/i r^diedii j-^d\7d\2 r^di^dio r^disdg r^digds r~)di4rfl3 
+ °*-^d2<i3 ^didi, ^diidi ^diedn-^digd-r-^digde 

/^ 1 Of? r)tZi6dii T\disdi2 j-\digdi3 r^di^dio r^di4ds 

^di5 ^^"-^dada ^dids ^dgdj ^ dudi ^ di^dn 
pidndg y-1 1/in r^diedii T-,disdi2 j-idi^dio r^dndi3 

-^digdig ^rflS ^^'^^d2d3 ^dids ^diidi ^disde 

j-\di4dg p^digds f^ I OOR pidi^dn pidisdi2 T^di:^dio 

^dndig ^digdj '-'dll, + ^«>0^d2d3 ^d4d5 ^dtidi 

r\di4dg j^digds pidirdi3 ^ , r/^ r^diedn ptdigdi2 

^diedij-^disd-r -^digde ^ di^^ ^ 0\) ly ^^^^ ^dids 

r-jdi^dio r-,disdg r)di4di3 p.dirds /-i r/^ T^d\jd\g j-jdisdii 

^diidi ^diedn-^disde ^digdy'-^di^ 0\) ly ^^^^ ^d^d^ 

pidiQdg r~)di5di3 T-,digds T-jdi4di2 /^ 99/1 r^diedig 

^di4di5^dnd-r ^ dmdi ^ digda ^'^16 ^^^-^^3^4 

Kit Kldi DttdV Kfdl' KlZ CdJ) /211680 



